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Quantum Spin Hall Effect and Enhanced Magnetic Response by Spin-Orbit Coupling 
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We show that the spin Hall conductivity in insulators is related with a magnetic susceptibility 
representing the strength of the spin-orbit coupling. We use this relationship as a guiding principle 
to search real materials showing quantum spin Hall effect. As a result, we theoretically predict that 
bismuth will show the quantum spin Hall effect, both by calculating the helical edge states, and by 
showing the non-triviality of the Z2 topological number, and propose possible experiments. 
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Spin Hall effect (SHE) Q, H Q has been attract- 
ing much attention recently, partly due to potential use 
for semiconductor spintronics. Its remarkable feature is 
to induce a spin current without breaking time-reversal 
symmetry. One of the interestingproposals is the quan- 
tum spin Hall (QSH) phase 0, H II II S 13 , which is 
a 2D insulator with helical edge states. The edge states 
form Kramers pairs, with spin currents flowing oppositely 
for opposite directions of spins. The QSH phase can be 
regarded as a novel phase constrained by the Z2 topo- 
logical number / |(J . It is equal to a number of Kramers 
pairs of helical edge states modulo two. The QSH phase 
has I = odd, while the spin-Hall-insulator (SHI) phase 
|lCj . topologically equivalent to a simple insulator, has 
/ = even. It is surprising that insulators without order- 
ing, usually considered as featureless and uninteresting, 
can have a nontrivial topological QSH phase. Its non- 
triviality reveals itself e.g. in a critical exponent 
For its interest akin to the quantum Hall systems, an 
experimental observation of the QSH phase is called for. 

To search for candidates for the QSH phase among a 
vast number of nonmagnetic insulators, we need a guid- 
ing principle. In this paper, we propose that the mag- 
netic susceptibility would be a good measure, and we 
pick up bismuth as a candidate due to its strong diamag- 
netism. There are also other supporting clues: large spin 
splitting in the surface states of the 3D system, similar 
to edge states for the 2D systems, the crystal structure 
in the (lll)_plane similar the Kane-Mele model for the 
QSH phase Using a 2D tight- binding model, we show 
that this system has only one pair of edge modes. Its Zi 
topological number [f| is shown to be odd, i.e. nontrivial, 
which supports stability of the edge state. These aspects 
make bismuth a promising candidate for the QSH phase. 

To relate magnetic susceptibility with spin Hall con- 
ductivity (SHC), we derive here a spin-Hall analog of 
the Stfeda formula. The Stfeda formula 0, El tells us 
that in insulators the Hall conductivity a xy is expressed 
as a xy = f^^lfi, where N is the number of states below 
the chemical potential /i, and is the area of the system. 
A direct generalization to the SHC is to replace N by the 
spin s z , i.e. a s — AtOjU- This is physically reasonable 
from the following argument. Suppose we apply a mag- 



netic field in some region, linear in time. Then a change 
of the total spin inside the region is proportional to 
According to the Maxwell equation, the increasing mag- 
netic field induces a circulating electric field. Therefore, 
by interpreting the spin change as due to a spin Hall cur- 
rent by the electric field, the SHC is <j s = jjib U- We 
can justify this argument by explicit calculations. 

A key step is to use a "conserved" spin current . In 
most papers on the SHC, the spin current is convention- 
ally defined as J s — ^{v, s z }. However, in the presence of 
the spin-orbit coupling, the spin is no longer conserved: 
T = s z ^ 0. Namely, the rhs of the equation of continu- 
ity, dts z + V • J s = T, is nonzero, and the "conventional" 
spin current J B is not directly related with spin accumu- 
lation. Instead, Shi et al. defined a conserved spin 
current J s as follows. If the system satisfies J dVT = 0, 
as it does in a uniform electric field, one can write T as 
T = — V • P T . Thus a conserved spin current defined as 
J s = J s + P T satisfies d t s, + V • J s = 0, and is calculated 
for several models |l41 Il5|. 

To calculate the SHC for the conserved spin current, 
we consider an electric field with wavenumber q, and 
take the limit q — » ^ij. When we calculate a spin 
current flowing to the x-direction in response to an elec- 
tric field to the y-direction, we take the vector potential 
A = A„ 



c y, and the response is calculated as 
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\{[H : s z ], e lqx }, the first and the second term correspond 
to the conventional and spin-torque terms (<r^° and aZ^ 
in Eq. (10) of Ref. 01)' respectively. By a calculation 
similar to [lij. we get — af^ 
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where G± = (s — H ± 177) 1 and L z = cci^ — yv x is an or- 
bital angular momentum. The term <j s is proportional 
to ^(G+ - G?_). By noting G+ - G_ = -2vr^( £ - if), 

only the states at the Fermi energy contribute to a s . 
In insulators vanishes identically. On the other 

hand, the second term a^ 11 ^ is expressed as 
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Equation agrees with the above-mentioned physically 
expected form. This result ||3J) can be also written as 
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where g is the electron-spin ^-factor, M or b. is an orbital 
magnetization, and /is is the Bohr magneton. These for- 
mulae are a spin analog of the Stfeda formula [lflj . We 
note that a Stfeda formula for the SHC with the conven- 
tional spin current J s has extra terms involving s z , in 
addition to These terms arise from spin nonconser- 
vation. Hence it is natural that they do not appear in 
as we used the conserved spin current [l7j . We remark 
that the definition of the spin current is still controversial. 
Since the spin is not conserved, there is no unique defini- 
tion of the spin current. Because there is no established 
way of directly measuring the spin current, one way is 
to consider instead a measurable quantities such as spin 
accumulation at edges. The spin accumulation depends 
crucially on boundary conditions, and the conserved spin 
current may correspond to smooth boundaries . This 
point requires further investigation. 

We calculate a s for the model on the honeycomb lattice 
proposed by Kane and Mele |fj|. This model shows the 
SHI and the QSH phases, depending on parameters Xr, 
X v and Xso ■ We numerically evaluate the SHC by Eq. (QJ 
using the formula of orbital magnetization |l8|, [l9| , 
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where is a sum over occupied bands. The result is 
shown in Fig. [J The calculated SHC is a s ~ in the SHI 
phase and a s ~ in the QSH phase, except for the 
vicinity of the phase boundary. The SHC in the QSH 
phase, <7 S ~ — e/(27r), is interpreted as a fundamental 
unit e/(47r) times two, the number of the edge states. The 
quantization is exact when s z is a good quantum number, 
i.e. \r — 0, where the system is a superposition of two 
quantum Hall systems with a xy — ±e 2 /ft. Remarkably, 
even when s z is no longer conserved, the SHC remains 
almost quantized. Deviation from the quantized value 
is more prominent near the phase boundaries, which is 
attributed to smallness of the band gap. 

Therefore, materials with large susceptibility would be 
a good candidate for the QSH phase. Namely, if the sus- 
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FIG. 1: Spin Hall conductivity a a for the Kane-Mele model 
on the honeycomb lattice for various values of Xr and X v in 
the unit of t. Aso is fixed as Q.QQt. As it is symmetric with 
respect to A_r — > — \r, we show only the result for negative 
Xr. The phase diagram is shown in the inset. Except for the 
vicinity of the phase boundary, o s ~ in the SHI phase, and 
<j s ~ —e/(2n) in the QSH phase. 



ceptibility is large, <j s should be large. FigureQ]then sug- 
gests that the system should be either in the QSH phase 
or near the phase boundary to the QSH phase. From this 
reason, we pick up some semimetals and related materi- 
als with large diamagnetic susceptibility, among which 
are bismuth and graphite. 

Bi crystal has a rhombohedral structure, with trigo- 
nal symmetry around the (111) axis. Bi is a semimetal 
with a small hole pocket at the T point, and three elec- 
tron pockets at the L points. Its strong diamagnetism 
has been studied experimentally and theoretically. It is 
theoretically attributed to massive Dirac fermions at the 
L and T points [2Cj. These Dirac fermions gives a dia- 
magnetic susceptibility which is enhanced as logarithm of 
the small energy gap |2("| . Even when the Fermi energy 
is in the gap, this picture survives, and the susceptibility 
becomes even larger. Such Dirac fermions contribute to 
anomalous Hall effect and the SHE. Hence, it is no won- 
der such enhanced diamagnetic susceptibility implies an 
enhanced charge/spin Hall conductivity. 

Because the QSH phase is in 2D, we have to make Bi 
two-dimensional, such as thin films and quantum wells. 
Such confinement discretizes the perpendicular momen- 
tum, and tends to open the gap. It was theoretically 
proposed that by making the Bi film thinner, it turns 
from semimetal to semiconductor [2lL l22j . Experiments 
show that the gap may open in thinner samples, whereas 
the gap is obscured by carrier unbalance between holes 
and electrons [23j. Remarkably, the lattice structure of 
Bi in the (111) plane resembles the Kane-Mele model ja]. 
The crystal can be viewed as a stacking of bilayers along 
the [111] direction. The inter-bilayer coupling is much 
smaller than the intra-bilayer one, and the LEED analy- 
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sis showed that the (111) surface of Bi is terminated with 
an intact bilayer [24| • Crystal structure of a single bilayer 
(Fig. Ufa)) consists of two triangular sublattices located 
in different layers. This honeycomb-like lattice structure 
is a key for a nontrivial Zi topological number; a crystal 
structure with high symmetry (e.g. square lattice) favors 
a trivial Z 2 topological number. 

We demonstrate that 2D singlc-bilayer bismuth has a 
pair of helical edge states carrying spin currents with 
opposite spins. Furthermore, we show that the Z 2 topo- 
logical number is odd. For these purposes, we use the 3D 
tight-binding model |25| which well reproduces the band 
structure, and truncate the model by retaining only the 
hoppings inside the bilayer. The resulting 16-band model 
is regarded as a multi-orbital version of the Kane-Mele 
model. We first calculate the band structure for a sin- 
gle bilayer (Fig. [5Ja)) for a strip geometry. The result 
is shown in Fig. [2[b), where projected bulk bands are 
shown in gray. The figure shows four edge states con- 
necting between the bulk conduction and valence bands. 
They correspond to one Kramers pair of edge states, sug- 
gesting nontrivial (odd) Z 2 topological number. The spin 
Chern number [2(j is calculated as C sc = —2, which is 
consistent with the existence of one pair of edge states. 
We also calculate the spin Hall conductivity from Eq. (0J 
and we get <r s ~ — 0.74 • This value is reduced from 
—2 • j- due to non-conservation of spin. 

To confirm that the 2D bismuth is in the QSH phase, 
we also calculate the Pfaffian Pf(fc) to calculate the Z 2 
topological number |(J. The bilayer system is inversion- 
symmetric, allowing Pf (fe) to be chosen real. The result is 
Fig.[2Ic), where Pf(fc) changes sign at the red curve, cor- 
responding to Fig. 2 (a) in [g; it implies odd Z 2 number. 
To further clarify the phase winding of Pf(fc), we break 
the inversion symmetry by adding small on-site energies 
±u, for the atoms on the upper and the lower layers, 
respectively. This may correspond to a heterostructure 
or a singlc-bilayer thin film on a substrate. The result 
is shown in Fig. |2td) for v — 0.2eV. There is only one 
vortex for the phase of Pf(fc) in the half BZ, which en- 
sures the odd Z2 number, corresponding to Fig. 2(b) in 
@. We note that the zeros of the Pfaffian do not fol- 
low the threefold rotational symmetry. It is because the 
Pfaffian is not covariant with respect to unitary trans- 
formation of the Hamiltonian, and depends on a choice 
of the unit cell. Because this QSH phase is protected 
by topology, it cannot be broken unless the valence and 
conduction bands touch at the same wavenumber and 
the direct gap closes. Therefore, even though the 2D 
tight-binding model might not reproduce quantitatively 
the real band structure, the nontriviality of the Z2 in- 
dex is more robust. This nontrivial Zonumber guaran- 
tees stability of the helical edge states |7J, |q • Only when 
the Z 2 topological number is odd, the edge states are 
stable against single-particle backscattering and (reason- 
ably weak) two-particle backscattering, whereas for even 




FIG. 2: Calculation on the bilayer tight-binding model of bis- 
muth in the (111) plane, (a) the crystal structure of the (111)- 
bilayer bismuth. The upper and lower layers are denoted by 
red and blue, respectively. The solid and broken lines repre- 
sent intralayer and interlayer hoppings, respectively, (b) the 
calculated band structure for the strip geometry with 20 sites 
wide. The gray region is the bulk bands, while the red curves 
are the states calculated for the strip. All states are doubly 
degenerate. Zeros of the Pfaffian Pf(fc) in the Brillouin zone 
are shown in red for (c) the inversion-symmetric (v = 0) and 
for (d) the inversion-asymmetric (v — 0.2eV) cases. 



topological number, the edge state will be gapped in gen- 
eral |7J,|3]. For example, in a bilayer antimony, we found 
that there are two pairs of edge modes, and the Z2 num- 
ber is even. Thus the edge modes in 2D antimony will 
be fragile against opening a gap. 

We now discuss a multilayer Bi thin film. Bismuth is 
suitable for pursuit of quantum size effects. A long mean- 
free path ( / ~ jim-mm) |27j , large electron mobility up 
to 10 6 cm 2 /Vs at 5K |2j|, and a small effective mass of 
electrons make the semimetal bismuth a good material 
to see quantum-size effects even at room temperature 
[2^|. Furthermore, Bi thin films can be synthesized with 
good quality. For example, a 10 /im-thick film shows 
magnetoresistance with a factor of few thousands at 5K 
and a factor of 2-3 at room temperature [28L Efflj l . 

By stacking N bilayers, the direct gap never closes, 
and each edge mode is topologically protected; the num- 
ber of pairs of edge modes becomes N. When the film 
becomes thicker, the helical edge states are expected to 
evolve to 2D surface states on a 3D bulk Bi, which has 
been studied by the angle-resolved photoemission spec- 
troscopy 0, llil HH ■ The surface states with large spin 
splitting is observed which manifests spin currents 
carried by such surface states. If the film thickness D is 
less than the mean free path I, e.g. D ~ 10 /im 
each of these edge modes has a quantized motion in the 
perpendicular direction. Hence, although the backscat- 
tering is relevant for even N, its effect on the edge states 
is almost negligible, and the system becomes gapless. By 
lowering the temperature or by increasing the disorder, 
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the system will eventually become the SHI or the QSH 
phases depending on whether N is even or odd. 

To observe the QSH phase experimentally, one way is 
to measure the spin current by an applied electric field. 
We note that the QSH phase does not show a quantiza- 
tion [(J, unlike the quantum Hall effect (QHE). On the 
other hand, surprisingly, the critical exponent v, govern- 
ing the localization length is found to be different be- 
tween the QSH and the SHI phases [Hi- Thus, the QSH 
phase can be established via measurement of v. In the 
QHE, v is determined experimentally by chan ging the 
magnetic field across the plateau transition 13.4 Isfij . 
because the change of the magnetic field controls the 
Fermi energy across the extended state. In the QSH 
phase, for example by a change of a gate voltage, one may 
be able to control the Fermi energy to the extended or 
localized states. To determine v, one has to see the range 
of the gate voltage A V„ showing nonzero a^ x by varying 
the sample size 36] or the temperature [3jj . From their 
critical exponents, v is calculated. 

Another way to establish the QSH phase is to 
observe the edge states by scanning tunneling mi- 
croscopy /spectroscopy, as has been used for graphite 
[371 138| . There is one important difference between the 
edge states of the graphite and the QSH phase. In 
graphene the existence of edge states crucially depends 
on the edge shape; the zigzag edge has edge states while 
the armchair edge does not. For a rough edge with por- 
tions of zigzag and armchair edges, only at the zigzag 
edges can the signal of edge states be seen [3^. In con- 
trast, the edge states in the QSH phase carry helical spin 
currents, ad circulate along the whole edge around, irre- 
spective of the details (e.g. the shape) of the edge. 

Besides bismuth, graphite is another material with 
anomalously large diamagnetic susceptibility |^. The 
spin-orbit coupling of graphite is small, and the dia- 
magnetism is mostly carried by orbital motion. In the 
same token as the SHE, the orbital -ang ular- momentum 
(OAM) Hall effect can be studied |4£j. The resulting 
OAM Hall conductivity is the susceptibility for the or- 



bital: 
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Since this involves only the 



n dB orb . 

orbital, the spin-orbit coupling is not required for it to 
be nonzero. Because the orbital susceptibility is largely 
enhanced in the graphite due to massless Dirac fermions, 
graphite will show large OAM current, 

In conclusion, we show that the SHC is directly related 
with a "spin-orbit" susceptibility which is a response of 
the orbital magnetization by the Zeeman field. We then 
propose that the magnetic susceptibility can be a good 
measure for search of quantum spin Hall systems. Wc 
theoretically predict that 2D bismuth will show the quan- 
tum spin Hall effect, because the number of pairs of he- 
lical edge states is odd and the Z 2 topological number is 
nontrivial. The sufrace states with large spin splitting in 
bulk bismuth, might be closely related with these edge 
modes. 
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